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This  paper  investigates  the  influence  of  electromagnetic  effects  on  the 
two  stream  instability  in  a  relativistic  electron  beam  propagating  through  a 
collisionless  plasma  channel.  The  analysis  is  carried  out  within  the 
framework  of  a  macroscopic  cold  fluid  model  in  which  beam  and  plasma  fluid 
element  is  in  a  laminar  flow.  Axisymmetric  electromagnetic  stability 
properties  are  calculated  for  the  case  in  which  the  equilibrium  beam  and 
plasma  density  profiles  are  rectangular.  Consistent  with  the  two  stream 
instability,  the  perturbed  fields  are  assumed  to  be  the  transverse  magnetic 
mode.  The  resulting  eigenvalue  equation  for  the  perturbed  axial  electric 
field  <5EZ  is  solved  to  give  a  closed  algebraic  dispersion  relation  for  the 
complex  eigenfrequency  a.  This  dispersion  relation  is  solved  and  it  is  shown 
that  the  electromagnetic  effects  have  a  strong  stabilizing  influence  for  a 
relativistic  electron  beam  with  y^  >  >  1,  where  Yb  is  the  relativistic  mass 
ratio  of  beam  electrons.  For  example,  the  critical  beam  current  for  instability 
is  proportional  to  the  electromagnetic  current  enhancement  factor  £  =  (y.  +  l)/2. 

y,  ~  D 

Thus,  the  critical  current  from  the  electromagnetic  calculation  increases  more 
drastically  with  the  beam  energy  than  that  from  the  electrostatic 
approximation.  *  -  ■  * 


I.  INTRODUCTION 


In  recent  years,  there  have  been  intense  research  activities  of  the 
equilibrium1’2  and  stability3"6  properties  of  relativistic  electron  beam 
propagating  through  a  background  plasma  channel  or  a  gaseous  medium.  Perhaps 
one  of  the  most  basic  instabilities  that  characterizes  a  relativistic  electron 
beam  propagating  through  a  collisionless  plasma  is  the  two  stream 
instability,6"9  which  results  from  the  relative  drift  motion  between  the  beam 
electrons  and  the  background  plasma  particles.  Although  the  two  stream 
instability  is  very  familiar  in  the  plasma  physics  community,  most  of  the 
previous  studies7"8  on  this  instability  have  been  limited  to  one  dimensional 
calculations.  In  recent  literature,6  Bogdankevich  and  Rukhadze  investigated 
the  two  stream  instability  of  a  relativistic  electron  beam,  including  the 
finite  radial  geometry  effects  on  stability  behavior.  Therefore,  they  were 
able  to  determine  the  limiting  beam  current  due  to  the  two  stream 
instability.  In  addition,  we  recently  developed  a  theory  of  the  two  stream 
instability  in  a  self-pinched  relativistic  electron  beam,9  finding  out  the 
critical  beam  current  for  instability.  However,  both  previous  calculations6,9 
have  been  based  on  the  electrostatic  approximation.  Although  this  is  a 
reasonable  approximation  for  a  mildly  relativistic  electron  beam,  we  expect 
significant  modifications  to  the  stability  behavior  when  the  beam  energy  is 
ultrarelativistic.  For  a  relativistic  beam,  the  electromagnetic  effects  on 
the  beam  dynamics  are  often  very  important.  In  this  regard,  in  this  paper,  we 
investigate  the  influence  of  electromagnetic  effects  on  the  two  stream 
instability  in  a  relativistic  electron  beam  propagating  through  a 
collisionless  plasma  channel. 
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The  analysis  in  this  article  is  carried  out  within  the  framework  of  a 
macroscopic  cold  fluid  model,  assuming  either  that  the  beam-plasma  fluids  are 
immersed  in  a  uniform  axial  magnetic  field  or  that  the  beam  is  self-pinched. 
Theoretical  model  and  the  equilibrium  properties  of  the  beam-plasma 
configuration  will  be  briefly  discussed  in  Sec.  II  for  rectangular  density 
profiles  of  the  beam  and  plasma  particles.  The  stability  analysis  presented 
in  Sec.  Ill  assume  axisymmetric  electromagnetic  perturbations  (3/38  =  0). 
Moreover,  since  the  unstable  mechanism  of  the  two  stream  instability  is  mostly 
due  to  fluctuations  of  the  axial  electron  field,  the  stability  analysis  in 
this  article  is  restricted  to  perturbations  of  the  transverse  magnetic  (TM) 
mode  polarization.  The  dispersion  relation  of  the  axisymmetric  transverse 
maynetic  mode  is  derived  in  Sec.  Ill  for  rectangular  density  profiles  of  the 
beam  and  plasma  particles. 

Assuming  long  axial  wavelength  perturbations  in  Sec.  IV,  we  apply  this 
dispersion  relation  to  two  specific  cases:  (a)  highly  magnetized  bounded 
plasma,  and  (b)  a  unmagnetized  self-pinched  electron  beam.  In  case  (a),  the 
dominant  two  stream  instability  results  from  the  relative  drift  motion  between 
the  beam  electrons  and  plasma  electrons.  One  of  the  most  important  features 
of  the  analysis  in  Sec.  IV  is  that  the  critical  beam  current  for  instability 
is  proportional  to  the  electromagnetic  current  enhancement  factor.  After  some 
analytical  algebraic  manipulation,  this  current  enhancement  factor  is 
approximately  given  by  [Eq.  (67)] 


for  the  equidensity  case  where  the  beam  and  plasma  electrons  have  the  same 

O 

number  density.  Here  Y^mc  is  the  characteristic  beam  electron  energy,  m  is 
the  electron  rest  mass  and  c  is  the  speed  of  light  in  vacuum.  In  the  limiting 
case  when  =  1 ,  this  enhancement  factor  recovers  the  electrostatic 
approximation  =  1.  Obviously,  the  critical  beam  current  for  instability 
increases  drastically  with  the  beam  energy  (y^),  clearly  demonstrating  that 
the  electromagnetic  effects  play  a  major  role  in  the  two  stream  stability 
behavior  for  a  relativistic  beam  with  Yb  >  >  1 .  In  Sec.  IV,  a  numerical 
investigation  of  the  dispersion  relation  is  also  carried  out  and  it  is  shown 
that  the  analytical  expression  of  the  current  enhancement  factor  in  Eq.  (67) 
underestimates  the  true  electromagnetic  current  enhancement.  Finally,  two 
stream  stability  properties  in  a  self-pinched  electron  beam  are  also  briefly 
discussed  in  Sec.  IV. 
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II.  THEORETICAL  MODEL  AND  EQUILIBRIUM  PROPERTIES 


1 1. A  Theoretical  Model 


In  the  present  analysis,  the  beam  and  plasma  particles  are  treated  as  a 
cold  fluid  immersed  in  a  uniform  axial  magnetic  field  Bq62.  In  the  limiting 
case  when  the  applied  axial  magnetic  field  vanishes,  the  beam  electrons  are  in 
a  self-pinched  equilibrium  where  the  self  magnetic  force  counterbalances  the 
centrifugal  force  resulted  from  the  beam  rotation.  Within  the  context  of  the 
macroscopic  cold  fluid  description,  the  equation  of  motion  and  the  continuity 
equation  for  the  beam  and  plasma  fluid  element  can  be  expressed  as 
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where  the  subscript  j  »  b,  i,  e  denotes  beams,  ions,  plasma  electrons, 
respectively,  P,  Is  the  mean  mechanical  momentum,  V  is  the  mean  velocity,  e, 
is  the  charge,  and  nj  is  the  density.  In  Eq.  (1),  E(x,t)  and  B(x,t)  are  the 
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electric  and  magnetic  fields  which  are  sel f-consi stently  determined  from  the 
Maxwell  equations,  that  is. 
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where  pgxt  and  Jgxt  are  externally  provided  charge  and  current  densities. 


I I. 8  General  Equilibrium  Properties 


An  equilibrium  analysis  of  Eqs.  (1)  -  (5)  is  carried  out  by  setting  3/9t  =  0. 


Thus,  the  equilibrium  properties  are  readily  obtained  from 
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where  n^(x),  V^(x),  P°(x),  E^(x)  and  8®(x)  are  the  macroscopic  equilibrium 
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quantities. 


We  introduce  cylindrical  polar  coordinates  (r,  9,  z)  with  z-axis 
corresponding  to  the  axis  of  symmetry;  r  is  the  radial  distance  from  the 
z-axis,  and  9  is  the  polar  angle  in  a  plane  perpendicular  to  the  z-axis.  For 
azimuthally  symmetric  particle  equilibrium  (9/59  =  0  and  3/9z  =  0)  characterized 
by  n^(r)  and  V  =  V^0(r)eQ  +  ViZ(r)ez>  ^  is  straightforward  to  show  from 
Eq.  (6)  that  the  functional  form  of  particle  density  profile  n^(r)  can  be 
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specified  arbitrarily.  Moreover,  the  deviation  from  equilibrium  charge  and 


current  neutrality  produces  a  radial  electric  and  azimuthal  magnetic  field 
that  influences  the  azimuthal  motion  of  particle  fluid.  It  follows  from 
Eq.  (7)  that  equilibrium  force  balance  in  the  radial  direction  can  De 
expressed  as 
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where  mj  for  beam  and  plasma  electrons  are  expressed  as  mb  =  ybm  and  mg  =  m, 
respectively,  and  the  equilibria  axial  velocity  profile  of  the  beam  electrons 
is  independent  of  r,  i.e.,  Vbz(r)  =  Sbc  =  const.,  which  relates  to  yb  by 
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yb  =  (l  -  3b)  .  We  also  assume  that  the  axial  velocity  of  ions  and  plasma 

electrons  is  zero,  although  this  restriction  can  be  easily  eliminated. 

However,  this  restriction  does  not  effect  on  the  electromagnetic  influence  of 
the  two  stream  instability.  In  Eq.  (10),  we  neylect  the  self  axial  magnetic 
field,  assuming  that  the  azimuthal  motion  of  particle  fluid  element  is 
non relativistic. 


With  oext  =  0,  the  equilibrium  radial  electric  field  is  determined  from 
Eq.  (9)  and  is  given  by 
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Similarly,  it  can  be  shown  that  the  equilibrium  magnetic  field  is  expressed  as 
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Substituting  Eqs.  (11)  and  (12)  into  Eq.  (10),  we  can  obtain  the  radial  force 
equation  explicitly. 


II.C  Sharp-Boundary  Equilibrium 


Although  the  formalism  outlined  in  Sec.  II. B  can  be  useful  to  investigate 
equilibrium  properties  for  a  broad  class  of  density  profiles,  for  purposes  of 
analytic  simplification  in  the  stability  analysis,  we  specialize  to  the  case 
of  a  sharp-boundary  equilibrium  in  which  the  particle  density  profiles  can  be 


expressed  as 
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In  Eq.  (13),  Rb  denotes  the  radius  of  the  beam  and  plasma  column.  As  a 
general  description,  we  also  assume  that  the  beam-plasma  system  is  bounded  by 
an  outer  conductor  with  radius  Rc.  However,  in  a  special  occasion  when 
Rc  *  » ,  the  system  becomes  unbounded.  For  convenience  in  the  subsequent 
analysis,  we  introduce 


fi  *  W  fe  =  W 


where  f^  and  fe  are  positive  constants, 


The  equilibrium  radial  electric  fie'ld  in  Eq.  (11)  is  expressed  as 
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Substituting  Eq.  (13)  into  Eq.  (12)  gives  the  azimuthal  maynetic  field 
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Laminar  rotational  frequencies  of  the  beam  and  plasma  fluid  elements  are 
determined  from  Eq.  (10)  by  substituting  Eqs.  (15)  and  (16)  into  Eq.  (10)  and 
carrying  out  a  straightforward  algebra.  Defining  the  laminar  rotational 
frequency  of  the  beam  electron  by  <ob  =  Vbe(r)/p*  we  ot)ta^n 
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where  wcb  =  eBQ/Ybmc  is  the  beam  electron  cyclotron  frequency,  u)pC)  =  4  ire  nb/ybm 

p  p  “1 

is  the  beam  plasma  frequency-squared  and  yp  =  (l  -  sb)  .  In  Eq.  (17),  the 
upper  sign  (u>b  =  id*)  corresponds  to  a  "fast1'  rotational  equilibrium,  and  the 
lower  sign  (u>b  =  <d“)  corresponds  to  a  "slow"  rotational  equilibrium. 
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Similarly,  we  also  obtain  the  rotational  frequency  of  the  plasma  ions  and 
electrons 
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respectively.  In  Eqs.  (18)  and  (19),  u>ci  and  uce  are  the  cyclotron 
frequencies  of  the  plasma  ions  and  electrons,  respectively,  and  n  =  y^n/m^ . 
From  Eqs.  (17)  -  (19),  we  require 
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for  radial  confinement  of  the  equilibrium.  The  inequalities  in  Eqs.  (20)  and 
(21)  assures  that  the  repulsive  space-charye  force  on  the  beam  and  plasma 
fluid  element  is  weaker  than  the  magnetic  focussing  force. 


III.  ELECTROMAGNETIC  STABILITY  ANALYSIS 


In  this  section,  we  linearize  Eqs.  (1)  -  (5)  assuming  electromagnetic 
perturbations  about  the  axi symmetric  equilibria  described  in  Sec.  II.  As 
indicated  in  the  introduction,  the  present  analysis  is  concentrated  on  the  two 
stream  instability  which  has  dominant  electric  field  polarization  along  the 
z-axis.  In  this  regard,  we  assume  that  all  perturbations  have  spatial 
dependence  only  on  the  z  coordinate,  according  to 

6<j>(x,t)  =  6i{/( r)exp{ i ( kz  -  wt)},  (22) 

where  6*(x,t)  represents  a  perturbed  quantity,  k  is  the  axial  wavenumber 
and  oj  is  the  complex  eigenfrequency.  The  assumption  of  the  axisymmetric 
perturbation  in  Eq.  (22)  simplifies  the  subsequent  stability  analysis 
considerably.  In  addition,  the  stability  analysis  in  this  article  is 
restricted  to  the  transverse  maynetic  mode  polarization  consistent  with  the 
electric  field  polarization  of  the  two  stream  instability.  Therefore,  the 
perturbed  field  equations  linearized  from  Eqs.  (3)  and  (4)  are  expressed  as 
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where  <5E  and  6B  are  the  perturbed  electric  and  magnetic  fields,  respectively. 


and  <5p(r)  and  6J(r)  are  perturbed  charge  and  current  densities,  which  are 


explicitly  defined  by 


Mr)  =  Iej5nj  . 
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5J(r)  =  XeJ(v96nJ  +  n9«V.), 
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respectively. 


The  linearized  equation  of  motion  and  the  linearized  continuity  equation 


for  the  beam  and  plasma  fluid  element  are  expressed  as 
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from  Eqs.  (1)  and  (2).  Before  proceeding  further  stability  analysis,  we 


restrict  the  subsequent  stability  study  to  the  following  cases:  (a)  in  the 


presence  of  an  applied  axial  magnetic  field  (Bq  *  0),  the  beam  density  is 


limited  to  satisfy 


8k  <  <  1. 
3  2b 
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where  v  is  Budker*s  parameter  of  the  beam  defined  by  v  =  e  N^/mc  ,  and  is 


the  number  of  beam  electrons  per  unit  axial  length,  or  (b)  in  case  Bg  =  0, 
self-pinched  electron  beam  propagating  through  a  pure  ion  channel  with  f.  <  1 
where  the  plasma  electrons  are  not  allowed  (f  =  0)  in  order  to  satisfy  the 
equilibrium  condition  in  Eq.  (21).  Within  the  context  of  this  restriction,  we 
neglect  the  terms  proportional  to  B0(r)  in  Eq.  (28),  substantially  simplifying 
the  subsequent  stability  analysis. 


After  a  straightforward  algebra  with  Eqs.  (28)  and  (29),  the  linearized 
fluid  calculation  gives 


■'(“  -  <1  4V  -  *  2“jJ  s5je  ■  if  ({ir  -  V?z«e/e).  l3» 

J 


•'(“  *  kV3z>  sije  *  +  2“j)  sijr  '  “• 

-1C-  - 


-  kvj*J  5V7lr(rnjiV)t,knj0s;jz*0 


where  the  relativistic  mass  ratio  y.  for  ions  and  plasma  electrons  is  unity, 
m.  =  Yhm  and  e  .  =  sgne  . . 
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From  Eqs.  (24)  and  (25),  it  is  straightforward  to  show  that  the 
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ination  6Ep  -  VjZ5B0/c  in  the  right-hand  side  of  Eq.  (31)  is  expressed  as 
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where  the  parameter  p  is  defined  by 
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After  a  tedious  algebraic  manipulation  which  makes  use  of  Eqs.  (31),  (32)  and 
(35),  we  obtain  the  coupled  equations 
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for  the  perturbed  radial  mean  fluid  velocities  <5V,r(r)  which  has  solutions  of 


the  form 


,  .  e.  (ca  -  kV?  )k  .  . 
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The  parameter  v?(<u,k)  in  Eq.  (40)  can  be  found  by  solving  Eqs.  (37)  -  (39) 

J 

simultaneously.  Note  from  Eq.  (40)  that  the  perturbed  mean  radial  fluid 

A  ^ 

velocity  5Vjr  is  expressed  only  in  terms  of  (d/dr)6Ez(r) . 

For  the  present  purposes  of  the  article  which  evaluates  the  influence  of 
electromagnetic  effects  on  the  two  stream  stability  behavior,  we  restrict 
present  stability  calculation  to  the  long  wavelength  perturbations 
characterized  by 


kRb  <  <  1. 


According  to  the  traditional  one  dimensional  two  stream  theory,  the  unstable 
k-values  in  a  typical  present  experiment  of  a  relativistic  electron  beam 
propagation  can  easily  satisfy  Eq.  (41).  In  the  context  of  Eq.  (41),  it  is 
not  necessary  to  explicitly  find  the  parameter  Vj(<i>,k)  in  Eq.  (40).  From 
Eq.  (33),  the  perturbed  mean  axial  fluid  velocity  is  also  given  by 


6V17  =  i  -J - - — r —  <SE  . 

J  'E (u  -  kvjZJ 


The  perturbation  in  density  in  Eq.  (34)  can  also  be  eliminated  in  favor  of 

A 

<5EZ .  Then,  after  a  straightforward  algebra,  it  is  shown  that  the  eigenvalue 
equation  in  Eq.  (23)  can  be  expressed  in  the  form 
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-  p Mi  -  1 


r]  <SE,  =  0, 


J  <“  -  «,*) 


2  2^ 

where  w  .  =  4ire  n./m,  is  the  plasma  frequency-squared  for  particles  of 

r  J  J  J 

9  9  9  9  9  9 

component  j  and  p  *  k  -  ur/c  is  defined  in  Eq.  (36).  Note  oj^  =  nf^^  dn 

2  2 
to  —  y  f  to  ,  . 

pe  b  e  pb 


Inside  the  beam-plasma  column  (0  <  r  <  Rb),  Eq.  (43)  can  be  expressed  in 
the  form 


-0. 


Outside  the  beam-plasma  column  (Rfa  <  r  c  Rc),  Eq.  (43)  reduces  to  Poisson 
equation  in  free  space,  that  is. 


7  Hr  (rHr  siz)  *  ’  u- 


2*  c  = 


R.  <  r  <  R  , 
b  c 


Choosing  a  right  solution  which  has  a  finite  value  at  r  =  0  and  properly 
matching  the  boundary  conditions  of  the  beam-plasma  surface  (r  =  Rb),  we 
eventually  obtain  the  dispersion  relation8’9 


u)V° 


[1  -  I  (1  -  -*&) 


C2k 


u,2.k2 

_£i_ 

27. 
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■J  tr 


J0<TV 

b  j^nv 


=  g(p) 


<0(PRc)  V*b>  '  K>Rb)  ^c* 

'  PKb  K0(pRc)  I0(pRb)  -  K0(pftb)  I0(pRc)  (47) 

where  Jq(x)  is  the  Bessel  function  of  the  first  kind  of  order  zero,  and  Iq(x) 

and  Kq{x)  are  modified  Bessel  functions  of  the  first  and  second  kinds, 

respectively,  of  order  zero,  and  the  geometrical  function  g(p)  represents  the 

right-hand  side  of  Eq.  (47).  The  "prime"  notation  in  Eq.  (47)  denotes 

derivatives  with  respect  to  the  complete  argument  of  the  Bessel  function, 

e.g.,  Iq ( x )  =  (d/dx)  IQ(x).  The  dispersion  relation  in  Eq.  (47)  is  one  of  the 

main  results  of  this  article,  which  can  be  used  to  investigate  the  two  stream 

stability  properties  for  a  broad  range  of  physical  parameters.  Particularly, 

the  influence  of  electromagnetic  effects  are  Incorporated  into  Eq.  (47) 

2  2  2  2 

through  the  parameter  p  =  k  -  <u  /c  .  Although  the  dispersion  relation  in 
Eq.  (47)  is  very  similar  to  that  of  the  electrostatic  approximation,8’9  the 
apparently  minor  modification  p  instead  of  k  makes  a  considerable  difference 
in  determining  the  critical  limiting  beam  current  due  to  the  two  stream 
instability,  as  will  be  seen  in  the  next  section. 


IV.  TWO  STREAM  STABILITY  PROPERTIES 


In  this  section,  making  use  of  the  dispersion  relation  in  Eg.  (47),  we 
investigate  two  stream  staDility  properties  for  long  wavelength  perturbations 
characterized  by  Eg.  (41).  The  dispersion  relation  in  Eq.  (47)  also 
represents  fluctuations  of  the  transverse  field  components  6Ep  and  6B0,  which 
may  drive  the  transverse  oscillation  of  the  beam  and  plasma  column.  For 
example,  it  is  more  involving  but  straightforward  to  obtain  the  dispersion 
relation  of  non-axi symmetric  electromagnetic  perturbations.  From  this 
dispersion  relation,  we  can  identify  typical  transverse  oscillations  such  as 
the  ion  resonance  instability10  resulted  from  the  transverse  oscillation  of 
the  beam  electrons  and  plasma  ions,  and  the  electron-electron  two-rotating 
stream  instability11  raised  from  the  relative  rotational  difference  between 
the  beam  electrons  and  plasma  electrons  which  is  unavoidably  introduced  by  the 
virtual  laminar  flow  assumption  in  the  cold  fluid  model.  The  instabilities 
associated  with  the  transverse  oscillations  are  other  important  issues  in  the 
electron  beam  propagation  experiments  and  they  must  be  separately  investigated 
with  more  attention.  For  the  present  purposes,  we  therefore  concentrate  on 
the  two  stream  instability,  which  is  a  characteristic  feature  of  fluctuations 
of  the  axial  electric  field.  However,  long  wavelength  stability  analysis 
somehow  eliminates  interference  of  the  transverse  oscillation,  thereby  leaviny 
behind  dispersion  relation  of  the  two  stream  instability,  which  will  be  used 
in  this  section  for  further  study.  For  the  long  wavelength  perturbations 
satisfying  Eq.  (41),  the  left-hand  side  of  Eq.  (47)  can  be  approximated  by 


uj2./y2 

p-Rr  > 


,2.}2  \ 

<b  4 
J  r 


=  2g(p) , 


- 


(48) 
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where  use  nas  Deen  made  of  |T|  <  1.  Within  the  context  of  the  assumption 

p2r2  <  <  j T j 2Rb2  <  1 ,  we  can  neglect  the  first  term  in  the  left-hand  side  of 
Eq.  (48).  However,  for  clarity  in  the  subsequent  stability  discussion,  we 
keep  this  term  without  effecting  final  results  in  this  article.  Before 
proceeding  stability  analysis  of  Eq.  (47),  we  investigate  general  properties 
of  the  two  stream  dispersion  function. 


IV. A  Analysis  of  Two  Stream  Dispersion  Function 


Here,  we  briefly  investigate  properties  of  a  typical  dispersion  function 
of  the  two  stream  instability,  which  is  expressed  as 


h2 

a  .  d 


u  ^  ^  -  i 


where  a,  b  and  c  are  positive  constants,  and  x  represents  complex 
eigenfrequency .  In  the  (f,x)  parameter  space,  f(x)  nas  its  local  minimum 


x  =  x>  =  5  (S^TT 


The  condition  for  instability  is  obtained  from  f(xp)  >  1  and  is  yiven  by 


C  <  b  [1  +  (a/b)  ^/-iJ 


?  2  .  . 

In  the  limiting  case  when  b  >  >  a  ,  the  minimum  point  xp  occurs  at 
Xp-  c .  After  a  careful  examination  of  the  dispersion  relation,  we  can  also 
show  that  the  maximum  growth  rate  occurs  at  the  axial  wavelength  satisfying 


K 


with  its  corresponding  solution 


/T  ra 


x  -  b  [  1  ±  i  — ^  ('25') 
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for  b4,  >  >  a  .  On  tne  other  hand,  in  the  limiting  case  b  <  <  a  ,  the 

2  1/3 

solution  x  -  (1/2)  (1  +  i  /T  )  (ab^/2)  of  Eq.  (49)  for  the  maximum  growth 
rate  occurs  at  the  wavenumber  z,  -  a. 


IV. B  Two  Stream  Instability  in  a  Magnetized  Beam-Plasma  System 


Two  stream  stability  properties  of  a  relativistic  electron  beam 
propagating  through  a  magnetized  plasma  channel  are  investigated  in  this 
subsection,  assuming  that  the  beam-plasma  system  is  bounded  by  an  outer 
conductor  with  radius  Rc.  Since  the  plasma  in  the  channel  is  strongly 
magnetized  (Bg  *  0) ,  not  only  ions  but  also  high  density  plasma  electrons  are 
allowed  in  the  channel,  still  satisfying  the  equilibrium  condition  in 
Eq.  (21).  Further  assuming  that  the  axial  wavelength  of  the  perturbation  is 
sufficiently  long  that 


P2R2c  <  <  1, 


it  is  straightforward  to  show  that  the  right-hand  side  of  Eq.  (47)  can  be 
approximated  by 


assuming  that  the  conducting  radius  Rc  is  reasonably  larger  than  the  beam 
radius  Rb.  Otherwise,  the  approximation  in  Eq.  (55)  is  failing.  Therefore, 
the  dispersion  relation  in  Eq.  (48)  is  approximately  written  by 


2 _ _  +  fk2  _  R  2 

CTO  -2  Kb 


0  2  0  “nk/Yw  (nf_.  +  YKfa)w^. 

=  (k2  “)r2  ( — 1  , 6  r  Pbj, 


(“  -  kShc)‘ 


which  can  be  numerically  solved  for  the  complex  eigenfrequency  u>. 


In  order  to  analytically  track  stability  properties  of  Eq.  (56),  we 


define 


k l  •  2/R l  m(Rc/Rb). 


Substituting  Eq.  (57)  into  Eq.  (56)  and  rearranging  the  terms,  we  can  express 
Eq.  (56)  as 

2  <A/y?  yKf  *2h 

M  pb  b  2,^4-Ebj  ■  1,  (58) 

k0  -  kebC^  “ 

2  2  2 

where  kg  kJ+  P  »  the  term  proportional  to  nf^  is  simply  dropped  since 

2  2  9  2  2 

n f i  <  <  Y^f g  in  a  typical  experiment.  Identifying  a  =  u)pbp“/YDky , 

2  2  2  2 

b  =  YbfeUpbp  /kg,  c  =  kfS^c  and  x  =  u,  Eq.  (58)  is  identical  in  structure  to 
Eq.  (49). 


From  Eq.  (51),  the  axial  wavenumber  k  of  the  instability  borderline  is 
approximately  given  by 


(whc)2  -  [l  +  — tttI  , 


in  which  the  function  f(x)  [defined  in  Eq.  (49)]  has  a  local  minimum  value  of 
unity  at 


u>  =  k6bc 


Ve'  +  1 


In  obtaining  Eqs.  (59)  and  (60),  we  have  assumed  that  p^  is  independent  of  the 
eiyenfrequency  w.  However,  it  is  a  function  of  to  in  reality.  Therefore, 
these  calculations  are  approximations.  For  better  understanding,  Eq.  (56) 
must  be  numerically  solved.  Substituting  Eq.  (60)  into  Eq.  (36),  it  is 
straightforward  to  show 


P2  -  k2/cth. 


at  the  axial  wavenumber  k  satisfying  the  stability  boundary  in  Eq.  (59).  In 
Eq.  (61),  the  electromagnetic  enhancement  factor  Ctn  is  defined  by 

(Ve’/3  *  'I' 

5th  ’  (62) 

Strictly  speaking,  the  theoretically  obtained  enhancement  factor  £  in 
Eq.  (62)  is  an  approximation.  The  meaning  of  "enhancement"  will  be  apparent 
in  the  following  discussion.  However,  we  emphasize  that  in  the  electrostatic 
approximation,  the  factor  r  in  Eq.  (62)  is  unity  for  arbitrary  value  of 


Yb.  Any  value  of  the  enhancement  factor  which  is  laryer  tnan  unity  is  the 
outcome  of  electromagnetic  effects. 


2  2 

Replacing  p  in  Eq.  (59)  by  k  and  rearranging  terms,  we  have  the 


relation 


(k8bc)2  =  Ybfeo>2b  [1  + 


Ve 


T77 


}  - 
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(63) 


for  the  stability  criterion.  In  obtaining  Eq.  (63),  use  has  been  made  of 
=  p  +  kj.  It  is  apparent  from  Eq.  (63)  that  the  electromagnetic  influence 
represented  by  increases  the  effective  transverse  wavenumber  from  kx  to 
k^^2  .  In  other  words,  rewriting  the  second  term  in  the  right-hand  side  of 
Eq.  (63)  by 


2.  .2  2  2Sbc2/t"  (Rc/Rb 

ith°C  R^/5 


V, 

we  observe  that  the  effective  beam  radius  reduces  from  Rb  to  Rb/Cth  ^  > 
thereby  enhancing  the  stabilizing  influence8’9  of  finite  radial  geometry 
effects  on  stability  behavior.  It  is  also  evident  from  Eq.  (63)  that  for 
instability,  the  beam  plasma  frequency  ojpb  must  satisfy 


25tb6bc2/Rb‘n(Rc/Rb) 


2  . 
(D  .  > 

pb 


Vet'  *  "Ve'”) 


7 


(64) 


which  is  necessary  condition  for  instability.  For  fixed  values  of  nb»Yb  and  f 
the  instability  condition  in  Eq.  (64)  cannot  be  satisfied  if  Rb  is 
sufficiently  small,  that  is,  finite  radial  geometry  effects  have  a  stabilizing 
influence  on  the  two  stream  instability. 


In  terms  of  the  magnitude  of  the  beam  current  I&  =  nbe8bcirRb,  it  is  straight¬ 
forward  to  show  that  the  instability  condition  in  Eq.  (64)  can  be  expressed  in 


the  equivalent  form. 


Jb  >  !crit’ 


where  the  critical  current  Icrit  for  instability  is  given  by 


A 

‘»(W 


:tn3b/fe 


In  Eq.  (66)  3  m^^/e  -  17000  amperes  is  the  Alfven  critical  current.  The 


electron  oeam  with  current  below  the  critical  current  is  stable  whereas  the 


beam  with  current  above  the  critical  current  is  unstable.  The  critical 


current  in  Eq.  (66)  is  proportional  to  the  current  enhancement  factor  £tb 


which  is  defined  in  Eq.  (62).  In  the  limiting  case  when  y  -  1 ,  the 


enhancement  factor  £  in  Eq.  (62)  recovers  the  electrostatic 


approximation  £th  *  1  for  arbitrary  value  of  fg.  In  order  to  dramatically 
demonstrate  dependence  of  £  on  the  beam  energy  yb,  we  simplify  Eq.  (62)  by 


for  the  equidensity  case  characteri zed  by  fe  =  1.  Obviously  from  Eqs.  (66) 


and  (67),  the  critical  beam  current  for  instability  increases  drastically  with 


the  beam  energy  yb. 


In  order  to  complete  the  stability  analysis  in  this  subsection,  we 


numerically  solve  the  dispersion  relation  in  Eq.  (56).  For  references  in  the 


ync< 

VAj 


subsequent  stability  analysis,  we  define  normalized  values  of  tne  critical 


beam  plasma  frequency-squared 


se  -  (-4^) 

c  crit 


2Bb/ltn(RC/RD) 
f=C  *  l/YRf  1/3)' 


for  the  electrostatic  approximation. 


2ttne^/fnRc/aD) 

sm  = - T7T  •  169 

V.O  ♦  l/Ve  ) 

from  Eq.  (64)  for  the  approximated  electromagnetic  calculation.  Note  that  an 
electron  beam  with  plasma  frequency  below  the  critical  plasma  frequency  is 
stable  whereas  tne  beam  witn  frequency  above  the  critical  plasma  frequency  is 
unstable.  We  also  remind  the  reader  that  the  critical  frequency  in  Eq.  (69) 
is  also  an  approximation.  From  a  numerical  calculation  of  Eq.  (56),  we  have 
obtained  correct  value  of  normalized  critical  plasma  frequency  -  squared  S 
which  represents  the  stability  boundary  in  (“pt,^/0  parameter  space. 

Shown  in  Fig.  1  is  plot  of  S  versus  Yb  obtained  from  Eq.  (5b)  for  Rc/Rb  =7.4 
and  fe  =  1 .  As  references,  we  also  present  Se  and  Sm  in  Fig.  1. 


The  current  enhancement  factor  in  Eq.  (62)  has  been  derived  from  the 
approximation  in  Eq.  (60).  Note  =  Sm/ Se  from  Eqs.  (68)  and  (69). 
Similarly,  we  numerically  determine  the  electromagnetic  current  enhancement 
factor  E,  from 


€  3  5/Se, 


»*? 


*8 


where  no  approximation  has  been  made.  Figure  2  presents  plot  of  C  versus 
obtained  from  Eq.  (56)  for  parameters  identical  to  Fig.  1.  The  theoretically 


obtained  is  also  presented  in  Fig.  2  as  a  reference.  Obviously  from 
Fig.  2,  we  note  that  the  theoretical  enhancement  factor  Cth  in  Eq.  (62)  is 
underestimating  the  true  electromagnetic  current  enhancement.  In  this  regard, 
the  true  critical  current  Icrit  for  instability  is  given  by  Eq.  (66)  where  the 
theoretical  value  of  is  replaced  by  the  numerical  evaluation  of  £.  As  a 
comparison  with  Fig.  2,  we  also  present  plot  of  C  versus  Yb  in  Fig.  3  for 
fe  =  10  and  Rc/Rb  =7.4  where  the  plasma  electron  density  is  ten  times  higher 
than  that  in  Fig.  2.  As  expected  from  Eq.  (62),  the  electromagnetic  current 
enhancement  increases  with  the  plasma  density. 


IV.C  Two  Stream  Stability  Properties  in  a 
Self- Pinched  Electron  Seam  ~~ 


As  a  second  example,  we  investigate  two  stream  stability  properties  in  a 
self-pinched  relativistic  electron  beam  propagating  through  a  pure  ion 
channel.  Because  of  Bg  =  0,  there  is  no  externally  provided  radial 
confinement  force  for  plasma  electrons,  thereby  not  allowing  the  plasma 
electrons  in  the  channel.  Moreover,  the  ion  density  in  the  channel  must 
satisfy  fj  <  1  in  order  for  a  radially  confined  equilibrium  for  ions  [see 
Eqs.  (18)  and  (19)].  We,  therefore,  concentrate  on  the  two  stream  analysis 
resulted  from  the  relative  drift  motion  between  the  beam  electrons  and  channel 
ions.  In  addition,  we  also  assume  that  the  beam  channel  system  is  unbounded 
with  Rc  ■*•  « .  Taylor  expanding  the  right-hand  side  of  Eq.  (47)  for  long 
wavelength  perturbations  [Eqs.  (41)]  and  making  use  of  Eq.  (48),  the 
dispersion  relation  can  be  expressed  in  the  approximate  form 
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are  stable,  the  beam  ion  two  stream  instability  is  stable  in  a  practical 
sense.9  Substituting  Eq.  (74)  into  Eq.  (73)  and  paralleling  the  derivation  of 
Eq.  (66),  we  obtain  the  critical  beam  current 


V.  CONCLUSIONS 


In  this  paper,  we  have  investigated  the  influence  of  electromagnetic 
effects  on  the  two  stream  instability  in  a  relativistic  electron  beam 
propagatiny  through  a  collisionless  plasma  channel.  The  equilibrium  and 
stability  analysis  (Secs.  1 1- IV)  was  carried  out  within  the  framework  of  a 
macroscopic  cold  fluid  model  in  which  the  beam  and  plasma  fluid  element  is  in 
a  laminar  flow.  Moreover,  axi symmetric  electromagnetic  stability  properties 
were  calculated  for  the  case  in  which  the  equilibrium  beam  and  plasma  density 
profiles  are  rectangular.  Consistent  with  the  two  stream  instability,  the 
perturbations  are  polarized  with  the  transverse  inaynetic  modes.  A  yeneral 
dispersion  relation  of  the  transverse  magnetic  modes  has  been  obtained  in 
Sec.  III.  This  dispersion  relation  has  been  applied  to  two  specific  cases. 

One  of  the  strikiny  features  of  the  stability  analysis  is  that  the  critical 
beam  current  for  instability  is  proportional  to  the  electromagnetic  current 
enhancement  factor  5,  which  increases  drastically  with  the  beam  energy  (yb). 
For  example,  the  equi density  case,  the  current  enhancement  factor  5  is 
approximately  given  by  5  =  (yfa  +  l)/2  for  the  beam  electron  -  plasma  electron 
two  stream  instability.  Thus,  for  a  relativistic  electron  beam  with 
Yb  >  >  1,  the  critical  beam  current  obtained  from  the  electromagnetic 
calculation  can  be  several  times  larger  than  that  from  the  electrostatic 
approximation.  We,  therefore,  conclude  that  the  influence  of  electromagnetic 
effects  on  the  two  stream  stability  behavior  plays  a  major  role  in  determining 
various  physical  parameters  for  a  long  distance  propagation  of  relativistic 
electron  beams. 
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FIGURE  CAPTIONS 


0  0  0 

Stability  boundary  in  (u)p|>R|/c  ,Yb^  Parameter  space  for 
Rc/Rb  3  ^ .4  and  f  =  1 .  S  is  numerically  obtained  from 
Eq.  (56).  Se  and  are  approximations  obtained  from 
Eqs.  (68)  and  (69),  respectively. 

Plot  of  £  versus  y^  obtained  from  Eq.  (56)  for  the 
parameters  identical  to  Fiyure  1.  As  a  reference,  £ 
is  also  presented. 

Plot  of  £  versus  yfa  obtained  from  Eq.  (56)  for  Rc/Rb  « 
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